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The Wave-making Resistance of Ships: a Theoretical and 

Practical Analysis . 
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§ 1. Introduction and Summary . 

The theoretical investigation of the total resistance to the forward motion 
of a ship is usually simplified by regarding it as the sum of certain 
independent terms such as the frictional, wave-making, and eddy-making 
resistances. The experimental study of frictional resistance leads to a 
formula of the type 

R,=/SV* (1) 

where S is the wetted surface, V the speed, / a frictional coefficient, and m 
an index whose value is about 1*83. 

After deducting from the total resistance the frictional part calculated from 
a suitable formula of this kind, the remainder is called the residuary resist¬ 
ance. Of this the wave-making resistance is the most important part; the 
present paper is limited to the study of wave-making resistance, and chiefly 
its variation with the speed of the ship. The hydrodynamical theory as it 
stands at present may be stated briefly. 

Simplify the problem first by having no diverging waves; that is, suppose 
the motion to be “ in two dimensions in space,” the crests and troughs being 
in infinite parallel lines at right angles to the direction of motion. Further, 
suppose that the motion was started at some remote period and has been 
maintained uniform. We know that, except very near to the travelling 
disturbance, the surface motion in the rear consists practically of simple 
periodic waves of length suitable to the velocity v of the disturbance. Let 
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Wave-making Resistance of Ships. 277 

cc be the amplitude of the waves, and w the weight of unit volume of water; 
then the mean energy of the wave motion per unit area of the water surface 
is \wa?. Imagine a fixed vertical plane in the rear of the disturbance ; the 
space in front of this plane is gaining energy at the rate \wo?v per unit time. 
But on account of the fluid motion, energy is supplied through the imaginary 
fixed plane to the space in front, and it can be shown that the rate of supply 
is \wa 2 u , where u is the group-velocity corresponding to the wave-velocity v. 
The nett rate of gain of energy is \wa 2 (v—u), and this represents the part of 
the power of the ship which is needed, at uniform velocity, to feed the 
procession of regular waves in its rear. An equivalent method of stating 
this argument is to regard the whole procession of regular waves from the 
beginning of the motion as a simple group; then the rear moves forward 
with velocity u while the head advances with velocity v, and the whole 
procession lengthens at the rate v—u. If we write B/y for the rate at which 
energy must be supplied by the ship, we call B the wave-making resistance, 
and we have 

B = \wa? (v—u)/v. (2) 

We notice that B is the wave-making resistance in uniform motion; it is 
only different from zero because u differs from v, that is, because the velocity 
of propagation depends upon the wave-length. 

In deep water, u is \v, so that B is \wa?. In the application of this to 
a ship at sea, it is assumed that the transverse waves have a certain average 
uniform breadth and height, and, further, that the diverging waves may be 
considered separately and as having crests of uniform height inclined at 
a certain angle to the line of motion; if the amplitude is taken to vary as 
the square of the velocity, it follows that B varies as v\ Several formulae 
of the type B = A'y 4 , or B = A^ 4 + B#, have been proposed ; although these 
may be of use practically by embodying the results of sets of experiments, 
they are not successful from a theoretical point of view. Becently <many 
such cases have been analysed graphically by Prof. Hovgaard ;* the general 
result is that a fair agreement may be made for lower velocities with an 
average experimental curve neglecting the humps and hollows due to the 
interference of bow and stern wave systems, but at higher velocities the 
experimental curve falls away very considerably from the empirical curve. 

The method used here consists in considering the ship, in regard to its 
wave-making properties, as equivalent to a transverse linear pressure 
distribution travelling uniformly over the surface of the water. Taking 
a simple form of diffused pressure system and making some necessary 

* W. Hovgaard, ‘ Inst. Nav. Arch. Trans., 5 vol. 50, p. 205, 1908. 
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278 Dr. T. H. Havelock. [Apr. 1, 

assumptions, we obtain an expression for the amplitude of the transverse 
waves thus originated, and for the resistance E, in which the velocity enters 
in the form er ^ v% ; this function is seen to have the general character of the 
experimental curves. Adding on a similar term for the waves diverging 
from bow and stern, and, finally, in the manner of W. Froude, an oscillating 
factor for the interference of these bow and stern waves, we find a formula 
for the wave-making resistance of the type 

E = ue~ l ^-\- /3 {1— 7 cos (m/v 2 )) e~ n ! v \ 

In this expression there are six adjustable constants; we proceed to reduce 
the number of these after transforming into units which utilise Froude’s law 
of comparison. We use the quantity c, defined as 

(speed in knots)/ (length of ship in feet), 

and we express the resistance in lbs. per ton displacement of the ship. An 
inspection of experimental curves, and other considerations suggest that the 
quantities l , m, n may be treated as universal constants ; with this assumption, 
a three-constant formula is obtained, viz,, 

E = ote- 2 ^ c * + /3 { 1-7 cos (10*2/c 2 )} e~ 2 ^, (3) 

where the constants a, /3, 7 depend upon the form of the ship. 

We then treat (3) as a semi-empirical formula of which the form has been 
suggested by the preceding theoretical considerations ; several experimental 
model curves are examined, and numerical calculations are given which show 
that these can be expressed very well by a formula of the above type. 

Since the constant a is found to be small compared with /3, it is not 
allowable to press too closely the theoretical interpretation of the first term, 
especially as the experimental curves include certain small elements in 
addition to wave-making resistance. If we limit the comparison to values 
of c from about 0*9 upwards, it is possible to fit the curves with an 
alternative formula of the type 

E = /3 {I -7 cos (10-2/c 2 )} 
and some examples of this are given. 

The effect of finite depth of water is considered, and a modification of the 
formula is obtained to express this effect as far as possible. Starting from 
an experimental curve for deep water, curves are drawn, from the formula, 
for the transverse wave resistance of the same model with different depths ; 
although certain simplifications have to be made, the curves show the 
character of the effect, and allow an estimate of the stage at which it becomes 
appreciable. 

In the last section the question of other types of pressure distribution is 
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1909.] The Wave-making Resistance of Ships . 279 

discussed, and one is given in illustration of the wave-making resistance of 
an entirely submerged vessel. 


§2. Pressure System travelling over Deep Water. 

It is known that a line pressure-disturbance travelling over the surface 
of water with uniform velocity v at right angles to its length gives rise to a 
regular wave-train in its rear of equal wave-velocity.* Take the axis of x 
in the direction of motion and let the pressure system be symmetrical with 
respect to the origin and given by p= f(x); suppose that f(x) vanishes 
for all but small values of x } for which it becomes infinite so that 


>00 

f(x ) dx = P. 


The regular part of the surface depression rj due to this 


integral pressure P practically concentrated on a line is given by 


2<yP . ax 
V = - JL o Sm S- 

WV J V z 


(4) 


The part of the surface effect which is neglected in this expression consists 
of a local disturbance symmetrical with respect to the origin and practically 
confined to its neighbourhood. 

If we suppose P constant, the amplitude in the regular wave-train and 
the consequent drain of energy due to its maintenance diminish with the 
velocity. 

To obtain results in any way comparable with practical conditions it is 
necessary to suppose the pressure system diffused over a strip which is not 
infinitely narrow. 

An illustration is afforded by taking 

y=/(a,)= I?i?’ < 6 > 

where a is small compared with the distances at which the regular surface 
effects are estimated. This type of pressure distribution is shown in fig. 1. 



* For a discussion of the wave pattern, see Lamb, 4 Hydrodynamics,’ § 241 et seq .; or 
Havelock, 4 Boy. Soc. Proc.,’ A, vol. 81, p. 398, 1908. 
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The effect of thus diffusing the pressure system is expressed by the 
introduction of a factor <£ (k) into the amplitude of the regular waves, where 
2tt j k is the wave-length and 





( 6 ) 


Using (5) in (6), we find 

<p (/c) = Ve~ aK = P^" a ^/ uS . 
Hence the amplitude of the waves is given by 


a = erwl*\ 


(7) 


Further, since tc = v 2 jg , the group velocity n = d(fcv)/dic = Hence 
the wave-making resistance E is given by 



( 8 ) 


We have to examine the variation of these quantities with the velocity^ 
under the supposition that the pressure system is due to the motion of a 
body either floating on the surface or wholly immersed in the water. The 
pressures concerned being the vertical components of the excess or defect 
due to the motion, it seems possible to assume as a first approximation that 
P varies as v 2 ; this is the case in the ordinary hydrodynamical theory of 
a solid in an infinite perfect fluid, and a similar assumption is also made 
in the theory of Froude’s law of comparison. This being assumed, we find 


a = A e~ a 9l v \ E = Be~ 2a ^ v \ 


(9) 


We see that both the amplitude and the resistance increase steadily from 
zero up to limiting values. 

If we draw the curve representing this relation between E and v, there is 
a point of inflection when 



or v 2 = ga . 


( 10 ) 


Writing v' for this velocity, we see that d^Rjdv increases as the velocity 
rises to v' and then falls off in value as the velocity is further increased. 

We can write the relation now in the form 


E = Be~WW. 


(ii) 


The character of this relation is shown by the curve in fig. 2, which 
represents the case 

E = 315e"f( 2 W (12) 

E being in tons, and Y in knots. 
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1909.] The Wave-mahing Resistance of Ships. 281 

The values of the constants in (12) have been chosen for comparison with 
an experimental curve of residuary resistance given by E. E. Eroude ;* it 
was obtained from model experiments and by means of the law of 



corresponding speeds and dimensions the results were given for a ship 
(model A) of 4090 tons displacement and 400 feet length. The actual curve 
is given in fig. 4 and is discussed more fully later; we neglect for the 
present the undulations which are known to be due to the interference of 
the bow and stern wave systems, and we consider a fairly drawn mean 
experimental curve denoted by E'. Table I shows a comparison of the 
values of E' with those of E calculated from the formula (12). 


Table I. 


Y. 

R. 

R'. 

10 

0-02 

1-8 

14 

2 

4 

18 

14 

16 

22 

38 

39 *5 

26 

70 

70 

30 

106 

107 

34 

132 

130 

38 

167 

166 

42 

176 

176 

46 

196 

192 


From this comparison we see that the point of inflection given by V' 
orresponds to the point at which the slope of the mean experimental curve 
* R. E. Froude, ‘Inst. Nav. Arch. Trans.,’ vol. 22, p. 220, 1881. 
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282 Dr. T. H. Havelock. [Apr. 1, 

begins to fall off. This effect is general in residuary resistance curves; we 
see that it is really an interference effect, the character of the curve being 
due to the mutual interference of the wave-making elements of the pressure 
system. Superposed on the mean curve we have a further interference effect 
due to the combination of two systems, the bow and stern systems. 

From Table I we infer that the mean curve agrees well with the calculated 
values E from about 18 knots upwards, but at the lower speeds the values of 
E are much too small; this suggests the addition of a term to represent the 
effect of the diverging waves. 

§3. Diverging Wave System. 

In the example considered above, the calculated values of E are much too 
small at the lower velocities. This might have been expected; for we 
obtained (12) by the consideration of line-waves on the surface, that is waves 
with crests of uniform height along parallel infinite lines. But the model 
experiments correspond more to a point disturbance travelling over the 
surface, with the formation of diverging waves as well as transverse waves. 
In fact, W. Froude* infers from his experimental curves that the residuary 
resistance at the lower velocities is chiefly due to the diverging wave system, 
on account of the absence of undulations; for the latter signify interference 
of the transverse systems initiated by the bow and stern, and these become 
very important at the higher velocities. 

We have to add to (12) a term representing the diverging waves; the 
comparison in Table I suggests for this a term of the same type, e-F v// / v >', 
with V" much smaller than the corresponding velocity \ r/ for the transverse 
waves. With the data at our disposal we might then determine the various 
constants so as to obtain the closest fit possible; however, we can make the 
process appear less artificial by the following considerations. We know that 
the wave pattern produced by a travelling point source consists of a system 
of transverse waves and a system of diverging waves, the whole pattern 
being contained with two radial lines making angles of about 19° 28' with 
the direction of motion; a fuller investigation of the effects produced by 
a diffused source must be left over at present. In applying energy con¬ 
siderations as in the previous sections, the usual method is to suppose that 
the transverse waves form on the average a regular wave-train of uniform 
amplitude and uniform breadth; using the same approximation for the 
diverging waves we suppose that these form on the average a regular wave- 
train on each side, with the crests inclined at some angle 6 to the direction 

* W. Froude, ‘ Inst. Nar. Arch. Trans.,’ vol. 18, p. 86, 1877. 
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1909.] The Wave-making Resistance of Ships . 283 

of motion of the disturbance. Then the velocity of the diverging wave- 
trains normally to their crests is Ysin#. Now the same features of the 
ship are responsible for the character of both transverse and diverging 
waves; then if V' is the velocity at which there is a point of inflection in 
the resistance curve for the transverse waves, the suggestion is that Y' sin 9 
is the corresponding velocity for the diverging waves. Taking as a first 
approximation the angle given above, viz., 19° 28' or sin” 1 ^, we test now 
a formula of the type 

E = A4-WFW+Be-WW, (13) 

For the particular example already used (Froude, Ship A) we take V' 
equal to 26 knots, and determine A, B from two values of Y. We obtain 
thus 

E = 4*5e“t(26/3V)* + 297^”^ 26 W (14) 

With this formula we find as good an agreement as before at the higher 
velocities, and we have now at lower velocities the comparison in Table II:— 


Table II. 


y. 

U. 

K'. 

10 

1*6. 

1*8 

14 

4*1 

4 

18 

16 ‘5 

16 

22 

40 

39*5 


In calculating from (14) we find that the two terms both increase 
continually; at low velocities the second term is practically negligible, then 
at about 15 knots the two terms are of equal value, and after that the 
transverse wave term becomes all important. 

It must be remembered that the experimental curve was obtained from 
tank experiments, and it is possible that the width of the tank may have an 
effect on the relative values of the transverse and diverging waves. It 
would be of interest if experiments were possible with the same model 
in tanks of different widths; if the methods used in obtaining (14) form 
a legitimate approximation, the effect might be shown in the relative 
proportions of the two terms—provided always that one can make a suitable 
deduction first for the frictional resistance, and can then separate out the 
relatively small effects of the diverging waves, the eddy-making and other 
similar elements. 
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§ 4. Interference of Boiv and Stern Wave-trains . 

The cause of the undulations in the resistance curves was shown by 
W. Froude to be interference of the wave system produced by the bow (or 
entrance) with that arising at the stern (or run). His experiments on the 
effect of introducing a parallel middle body between entrance and run 
confirmed his theory, which may be stated briefly. Let the wave-making 
features of the bow produce transverse waves which would have at 
a breadth b an amplitude a ; owing to the spreading out of the transverse 
waves they will be equivalent to simple waves at the stern of smaller 
amplitude ha, at the same breadth b. Let a' be the amplitude there of the 
waves produced by the stern. Then in the rear of the ship we suppose there 
are simple waves of amplitude ha superposed upon others of equal wave¬ 
length of amplitude a'. At certain velocities the crests of the two systems 
coincide in position, giving rise to a hump on the resistance curve; and at 
intermediate velocities there are hollows on the curve owing to the crests of 
one system coinciding with the troughs of the other. 

In developing a form for the resistance, subsequent writers have generally 
taken E proportional to an expression of the form a 2 -\-a' 2 ~\-2haa' cos (mghjv 2 '), 
where L is the length of the ship. This means that the bow is supposed to 
initiate a system of waves with a first crest at a short distance behind the 
bow, and that similarly the stern waves have their first crest shortly after 
the stern; the length mL is the distance between these two crests, and is 
called the wave-making length of the ship. The determination of a value 
for m appears to be doubtful, but from interference effects it is said to vary 
for different ships between the values 1 and 1*2. 

It has seemed desirable here to follow more closely the point of view in 
W. Froude’s original paper already quoted.* We regard the entrance of the 
ship as forming transverse waves with their first crest shortly aft of the bow, 
and the run of the ship as forming waves with their first trough in the 
vicinity of the .middle of the run. It is suggested that this distance between 
first crest and first trough, in practice found to be about 0*9L, should be 
taken as the “wave-making distance”; the cosine term in the formula 
is then prefixed by a minus sign instead of a positive sign. We return to 
this point later; we first work out a definite simple illustration in “ two- 
dimensional waves,” and then build up a more complete formula for 
comparison with experiment. With the same notation as in § 1, let the 
pressure system be given by 


P 


-/M-U - P 2“ 2 1 

" 1 K ) ~ 7T U 2 +(® - ¥f +(® ■+ iO 2 J 


* W. Froude, loc. cit. ante , p. 83. 


(15) 
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1909.] The Wave-maJcing Resistance of Ships. 285 

This indicates two pressure systems, one of excess and the oth§r of defect 
of pressure; each distribution is of the type already used, and their centres 
are separated by a distance l. Fig. 3 shows the character of the disturbance. 



In the rear of the whole disturbance there is interference between the 
regular wave-trains due to the two parts. With the same methods as before 
we find that the resulting waves are given by 

v =Ml e - °*/°> S in gfr-** )— e~W sin l^+ M l 

WV 2 V 2 WV 2 V 2 

= ~p e ~ a9 ' vi {(Pi-P»)008^8in^-(P 1 +.P 8 ) sin J^cos^f}. (16) 

Hence the average energy per unit area is proportional to 
v -± e -2aglv> { Pl 2 + p 2 2_2P 1 P 2 COS (gl/v 2 )}. 

How, assuming as before that Pi and P 2 vary as v 2 , we find that as regards 
variation with the velocity the effective resistance K, which is the expression 
of the energy required to feed the wave-trains, is given in the form 

K = {A 2 + B 2 —-2AB cos (< glfv 2 )} <r 2a ^\ (17) 

A more general expression might have been obtained by taking two 
quantities and a 2 in (15), corresponding to some difference in wave-making 
properties of entrance and run; this would have led to different exponential 
factors being attached to the bow and stern waves. However, we find (17), 
with a common exponential factor, sufficiently adjustable for present 
purposes. 

In Froude’s experiments in 1877 the effect of inserting different lengths 
of parallel middle body between the same entrance and run was examined; 
it was found that a hump in the residuary resistance curve corresponded to 
a trough of the bow waves being in the vicinity of the middle of the run, 
and a hollow to a crest being in that position. 
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286 Dr. T. H. Havelock. [Apr. 1> 

For the model, Ship A, we have: Length = L = 400 feet; entrance = 
run = 80 feet. 

Hence, in this case we may take, in formula (17), / as approximately 
360 feet. We notice that this gives l = 0*9L; and in subsequent com¬ 
parisons, instead of leaving l to be adjusted to fit the experimental curve, 
we find there is sufficient agreement if we fix it beforehand as 0*9 of the 
length of the ship on the water-line. 

Compare, now, the length l with the ordinary “ wave-making length ” of 
the ship; the latter is written as mL and is defined as the distance between 
the first regular bow crest and the first regular stern crest. From the present 
point of view (17) gives 

mL = Z + JX or m = 0*9-h^X/L, (18) 

where X is the wave-length in feet of deep-sea waves of velocity v ft./sec. 

Calculating from this formula for Ship A, and writing Y for velocity in knots 
(6080 feet per hour), we obtain Table III. 

We see that the statement that m lies between 1 and about 1*2 would 
hold for this ship if it were measured for ordinary speeds between about 
14 and 22 knots. 


Table III. 


V. 


m. 

10 

65*5 

0*97 

14 

110 

1-03 

18 

180 

1-12 

22 

270 

1 -24 

26 

362 

1-35 

30 

500 

1*5 


We proceed now to modify (14) by introducing into the second term a 
factor I —7 cos(#7/^ 2 ). With l = 360, we find gljv 2 is approximately 
4080/Y 2 , with V in knots; further, from one value from the experimental 
curve we obtain 7 = 0T2. Thus for Ship A we have K in tons given by 

E = 4*5^t(26/3V)« + 2 97 {1-0T2 cos (4080/Y 2 )} rtWY (19) 

Table IY shows some calculated values for E, and these are represented in 
fig. 4 by dots; the continuous curve is the experimental residuary resistance 
curve given by Froude, that is, the total resistance less the calculated 
frictional part. 

It is the custom to give the results of model experiments in the form 
of a fair curve, so that the positions of actual readings and the possible 
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1909.] The Wave-making Resistance of Ships. 287 



10 15 20 V Knots 25, 30 35 40 45 50 


error are not known. The interrupted curve is a curve R = AY 4 sketched in 
for comparison. 


Table IY. 


Y. 

R. 

Y. 

R. 

10 

1*5 

30 

102 

14 

4*2 

34 

142 

18 

15 

38 

171 

22 

44 

42 

185 

26 

62 

46 

195 


§ 5. Comparison with Experimental Results. 

Before examining further model curves we must express the previous 
formula in a form more suitable for calculation; we use the system of units 
in which model results are now generally expressed. R is given in lbs. 
per ton displacement of the ship, while instead of the speed Y we use the 
ratio YIx/L, Y being in knots and L in feet; this is called the speed-length 
ratio, and we shall denote it by c. The advantage of these units is that they 
utilise Froude’s law of comparison; from the experimental curve between 
R and c we can write down at once the residuary resistance for a ship of any 
length and displacement at the corresponding velocity, provided the ship has 
the same lines and form as the model. Thus the constants which are left in 
the relation between R and e depend only upon the lines of the model, not 
upon its absolute size. At present we make no attempt to connect these 
constants with the form of the model, as expressed by the usual coefficients 
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288 Dr. T. H. Havelock. [Apr. 1, 

of fineness or the curve of sectional areas, or in other ways; we are concerned 
with the form of E as a function of c , and the constants are chosen in each 
case to make the best fit possible. 

First, as regards the exponential factor, we had g-^v'/vp with V' giving 
a point of inflection on the resistance curve; in the case of Ship A we had 
V' = 26, L = 400, so that c = 1*3. Now, it is just about this value of c 
that there is a falling off in most experimental curves, so that we try first 
c' = 1*3 for the point of inflection on the E, c curve. Then the exponential 
factor becomes er% c '/ c)2 , or 

■ Secondly, as regards the cosine term which gives the undulations, we had 
o,os(gl/v 2 ); we have decided to put l = 0*9L, so that we have 

fi = °% L /(S V ) 2 = approximately. 

Hence the previous relation for E reduces to the following general form: 

E = aer 2 * 53/9c2 + /3 (1 7 cos W2 / c 2 ) e~ 2 'W c \ (20) 

where E is in lbs. per ton displacement, and a, / 3 , y depend upon the form of 
the model. 

There are humps on the curve when 10*2c ~ 2 is an odd multiple of t r, 
hollows when it is an even multiple, and mean values when it is an odd 
multiple of For facilitating calculation, some of these positions are given 
in Table Y; and, for the same reason, values of the exponentials and the 
cosine factor are given in Table YI. 


Table Y. 


Humps 

Means ... 

— 

2*54 

1 *8 

1 *47 

— 

1 -13 

1-04 

0*96 

— 

0*85 

0*8 

0*76 


Hollows 

00 

— 

— 

— 

1 -27 

— 

— 

— 

0*9 

— 

— 

— 

0-73 


Values of c. 


Table VI. 


c. 

g—2‘53/9c 2 > 

e -2 -53/c*. 

cos (10*2/c 2 ). 

0-6 

0 

'460 

0 *0009 

+ 0-75 

0*8 

0 

•644 

0*019 

-0-97 

1 *0 

0 

•756 

0*080 

-0*71 

1 *2 

0 

•821 

0*172 

+ 0 -70 

1 *4 

0 

•866 

0-275 

+ 0*47 

1 *6 

0 

•896 

0-372 

-0*65 

1 *8 

0 

•916 

0-458 

-1 *0 

2 *0 

0 

•932 

0 -532 

— 0 *83 

2 -2 

0 

•943 

0-592 

— 0 *51 

2*4 

0 

•951 

0-644 

-0-20 

3 

0 

•970 

0*756 

+ 0*43 
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1909.] The Wave-making Resistance of Ships. 289 

We examine, now, some examples of experimental curves, comparing them 
with the formula (20) ; several of the curves and other data, in particular 
for II, III, and V, have been taken from the collection in Prof. Hovgaard’s 
paper already referred to, in which he essays to fit formulae involving V 4 or 
V 6 with the experimental curves. 

I. R. E. Froude , 1881, Ship A. 

Displacement = 4090 tons ; length = 400 feet; cylindrical 
coefficient = 0*694.. 

This is the case we have examined in the previous sections, so that we 
have only to change the numerical factors in (19) to cause E to be given in 
lbs. per ton displacement. We find the result is formula (20) with 

a = 2*46; /3 = 162*6; y = 0*12. 


II. W. Froude, 1877. 

Displacement — 3804 tons; length = 340 feet; cylindrical 
coefficient = 0*787. 

The last two data include the cylindrical middle body. The curve is 
given in fig. 5 ; it was constructed by Hovgaard from the data of Froude’s 
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290 Dr. T. H. Havelock. [Apr. 1, 

experiments, and these were such that it was possible to make a mean 
residuary resistance curve, the effects of bow and stern interference being 
eliminated. The curve is given as total residuary resistance in tons on a 
base of Y in knots. If we work in lbs. per ton, we find there is a very fair 
agreement with formula ( 20 ) if we take 

a = 2*24; /3 = 279*7; 7 = 0 . 

Probably a closer agreement could be obtained by further slight adjustment 
of a and /3. Pig. 5 shows a comparison of values of the total residuary 
resistance for the ship (in tons) ; the calculated values are indicated by. small 
circles. 

III. R W. Taylor , 1000 lbs. Model 
Length on water line = 20*51 feet; cyl. coeff. = 0*680. 

The experimental curve in this case is given as residuary resistance for 
the model in lbs. on a base of Y in knots. With the same notation as before 
we find 

a = 2 ; /3 = 136*6 ; 7 = 0*14. 

Putting these values in (20), we can calculate E in lbs. per ton, and hence 
Ei in lbs. for the model; fig. 6 shows the comparison between Ei and the 
corresponding values on the curve; the calculated values Ei are indicated by 
dots. 
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1909.] TTie Wave-making Resistance of Ships. 291 

IV. D. IF. rayZor, Model No. 892.* 

Displacement = 500 lbs.; length on water line = 20*512 feet; longitudinal 
coeff. = 0*68; midship section coeff. = 0*70. 

In this case the experimental curve is given as lbs. per ton displacement 
(R') on a base of speed-length ratio ( c ). In the same manner as before, fig. 7 
shows the comparison with the formula (20) when we take 

a = 2 ; 0 = 82*5 ; 7 = 0*14. 

Since the constant a is small compared with 0, one is not able to lay 
much stress on the meaning of the first term. For as the velocity functions 

50 

40 


30 


20 


IO 


•5 .7 C 9 1.1 13 1.5 1.7 1.8 

are of a suitable type, the constants possess considerable elasticity as regards 
fitting an experimental curve. For instance, if we omit values of c below 
about 0*9, it is possible to represent the previous curves fairly well by a 
formula 

R = 0 {1—7 cos (10*2/c 2 )} 

In the previous examples we took the value 1*3 for c'. In Case IV above 
we find now the values 

/3 = 87; 7 = 0*14; V = 1*3. 

For a similar curve taken from the same paper, viz., Model No. 891, dis¬ 
placement 1000 lbs., we find a good correspondence, except for slightly higher 
values near c = 1*1, with the values 

0 = 174 ; 7 = 0*14; c' = 1*4. 

* D. W. Taylor, Society of Naval Architects, New York, November 19 , 1908 . 

VOL. LXXXII.—A. 
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292 Dr. T. H. Havelock. [Apr. 1, 

V. I. I. Yates, Destroyer Model C* 

Displacement = 575 lbs.; length = 20 feet; cyl. coeff. = 0*529. 

The experimental curve is given in lbs. for the model on a base of Y in 
knots, and is a total resistance curve, that is, it includes the frictional 
resistance. The curve is reproduced in fig. 8. 



This curve is not analysed here so as to compare the residuary resistance 
with the formula (20), but it is included in order to draw attention to certain 
possible complications. It may be noticed that the curve is carried to a high 
value of the speed-length ratio c, and that it continues to rise more rapidly 
after about e = 2*3 than might be expected on the present theory. ISFow’ in 
the first place it is possible that the frictional resistance may account partly 
for this rise. The ordinary estimation of the frictional resistance assumes 
that it can be calculated separately from some expression like / SV r85 ; now 
the legitimacy of this is beyond doubt in all ordinary cases, but at high speeds 
it is possible that the form of the expression may change, or even that it 
may not be a fair simplification to divide the total resistance into simple 
additive components. 

In the second place a more important consideration must be taken into 
account, and that is the depth of the tank. Tor the experiments now under 
* I. I. Yates, Thesis, lQO 1 ?, Mass. Inst. Tech. U.S.A. See Hovgaard, loo. eit. ante. 
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consideration the depth of water in the tank is not known. The deepest 
experimental tank appears to be the U.S. Government tank at Washington, 
which has a maximum depth of about 14*7 feet. Now in that tank, with 
a 20-foot model, there would be a “ critical ” condition near the value 
c = 2*9 ; before and up to that point the residuary resistance curve would 
rise sharply and abnormally. This effect is discussed more fully in the next 
section, and curves are given in fig. 11, with which fig. 8 may be compared. 
It appears, then, as far as one is able to judge, that it is possible the 
resistance curve in fig. 8 is complicated by the effect of finite depth of the 
tank. 


6. The Effect of Shallow Water . 


We saw in the first section that the wave-making resistance R can be 
written in the form 


R = \wa 2 (v—u)/v, 


where u is the group-velocity corresponding to wave-velocity v. For deep 
water u = \v 9 and the formulae are comparatively simple. But for water of 
finite depth h the relation between u and v depends upon the wave-length 
(27r//e). We have 



u = (/cv) = (1 + 2/eA/sinh 2tch). 

a/c 


Consequently we find 



( 21 ) 


As v increases from zero to y/(gh), R diminishes from \wa? to 0, provided 
the amplitude remains constant. But as Prof. Lamb remarks,* the 
amplitude due to a disturbance of given character will also vary with the 
velocity. It is the variation of this factor that we have to examine in 
the manner used in the previous sections for deep water. 

If a symmetrical line-pressure system F(oj), suitable for Fourier analysis, 
is moving uniformly with velocity v over the surface of water, the surface 
disturbance 7) is given by 


ttwt] = £ ^ /cV<j> (k) sin /c {x + (y — Y) t} die 


0 Jo 


— J | dt J fcY(j)(/c) sin K {a:+(i? + V)£} d/c, (22) 


0 Jo 


0 Jo 



* H. Lamb, 4 Hydrodynamics, 5 p. 391, 1906. 


x 2 
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The method of evaluating these integrals approximately so as to give the 
regular wave-trains has been discussed in a previous paper and it is followed 
now in the case of finite depth.* We take, under certain limitations, the 
value of an integral such as 



to be the value of its principal group, viz., 



(22a) 


where tc 0 is such that g' (w 0 ) = 0. 

Now in the integrals in (22) we have to find successively two principal 
groups, first with regard to ic and then in the variable t ; and thus we may 
evaluate the amplitude factor in the resulting regular wave-trains. 

Tor water of depth h we may write 



The group with respect to k gives a term proportional to 

cos {tK?f (/c) + |-7r}, 
where k has the value given by 


/(«)+«/'(*)= -j- 


(23) 


From (22a), this introduces into the amplitude a factor 


l/v/p {2/ «+*/" (*)}]• 

Further, the group with respect to t occurs for 


(24) 


J {*#*/(«)}= 0 or /(«) = 0. 


Also we have in these circumstances 





* f+«f l (/+«/') 2 _ (VO 2 
t 2 2/ + V" t 2/ + a :/" t (2 If + «/0 ’ 


(25) 


t (2f' + */')• 


Hence from (22a), (24), and (25) the selection of the two groups adds to 
the amplitude a factor 1 j k/^k), where 


/(#) = 0=v- 


,IC 


* Havelock, ‘Roy. Soc. Proc.,’ A, vol. 81, p. 411, 1908. 
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Also if u is the group-velocity for wave-length 27 t/k and wave-velocity Y, 
we have, in this case, 


u = Jlc ^ V ' ) = Tk {^ - */(*)} = »-{/(*)+*/' («)}• 


Hence, since in the final value /(#) = 0, we have equal to v—u. Thus 

if k is the wave-length of the regular wave-trains in the rear of the 
disturbance, we find that they are given by 


, k: vd> (k) • 

7} = const, x . y - v -. sin tcx, 

v—u 


(26) 


where v = /y/^ tanh ^ ^1 + 

Hence for the amplitude a we have ^ 


2 kJi 


where 


sinh 2 tch} * 



Substituting now in (21) we obtain for the wave-making resistance, E 
proportional to 



If we take the same distribution of pressure in the travelling disturbance, 
namely, F (%) = Pa/7r (a 2 -i-% 2 ), we have (f> (/c) = Pe~ aK ; further, we may 
again assume that the pressure P varies as v 2 f so that we have the resistance 
in the form 



tanh kJi v 2 


with 


(27) 


tch gh ’ 


Considering E given as a function of v by these two equations, we see 
that E increases slowly at first and then rapidly up to a limiting value at 
the critical velocity after this point E is zero, for there is no value of 

k satisfying the second equation with v 2 jgh>l. 

Further, the limiting value of E at the critical velocity is finite, for we 
have 


1*5. 


« = o (1 — 2^/sinh 2 kK) 


We see that the E, v curve given by (27) is of the type sketched in fig. 9, 
We may compare this with some of the curves given by Scott Eussell for 
canal boats. The continuous curve in fig. 10 is an experimental curve of 
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296 Dr. T. H. Havelock. [Apr. 1, 

total resistance,* and the dotted curve is a parabolic curve inserted here to 
represent approximately the frictional resistance; the difference between the 
two curves represents the residuary resistance, and is clearly of the same 
type as the theoretical curve in fig. 9. 




We can obtain a better estimate of equation (27) by taking an experi¬ 
mental curve for a model in deep water, and then building up curves for 
different depths. We must first put (27) into a form suitable for com¬ 
parison with deep water results. 

Limiting the problem to one of transverse waves only, the formula (27) 
must reduce to R = Ae” 2 ’ 53 ^ 2 , for h infinite andc = (speed in knots)/^/(length 
in feet). 

Writing v' for v/^(gh) we find c 2 = ll‘3?/%/L; thus although the actual 
critical velocity does not depend upon the length of the ship but only on the 
depth of water, the speed-length ratio (e) has a critical value which is 
proportional to the square root of the ratio (depth of water)/(length of ship). 

In (27) we cannot fix any value of v or c and then calculate R directly; 
we must w^ork through the intermediate variable kJi. The equations may 
now be written as 

R = A (fch) 2 v%-P *7(1 - 2^/sinh 2 «h), (28) 

v'* = (tanh k1i)/k1i ; (3 ' = 0-218L/A; c 2 = lVSv' 2 h/L. 

With h infinite this reduces to the previous form for deep water with the 
same constant A, so that a direct comparison is possible. As the velocity v 
increases from 0 to \/(gh), & diminishes from oo to 0 ; we select certain 
values of kJi, calculate the values from tables of hyperbolic functions, and 
thus obtain the set of values in Table VII, writing m for 

(^)V 4 /(l-2^/sinh 2 Kh). 

* J. Scott Kussell, ‘ Edin. Phil. Trans./ vol. 14, p. 48, 1840. 
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1909.] The Wave-making Resistance of Ships. 297 

Table VII. 


Kh. 

v\J(gli). 

c 2 L/A 

m. 

— £ /cc? 2 . 

00 

0 

0 

1 -o 

2‘53 

10 

0*316 

1-13 

1 -o 

2-53 

6 

0*41 

1-87 

1-0 

2-53 

4 

0*5 

2-82 

1*005 

2*53 

2 

0*69 

5-42 

1-077 

2 -43 

1 

0-87 

8*57 

1 *287 

1-92 

O 

1-0 

11 -3 

1 -5 

0 


We consider now the experimental curve analysed in Case IV in the 
previous section, a model of 20*5 feet taken up to a value c = 1*8. Assuming 
that the influence of finite depth was inappreciable in this range, we have for 
deep water 

E = 2e- 2 ’ 53 / 9c2 + 82*5 {1 — 0*14 cos (10*2 /c 2 )} e ~^\ (29) 

We leave out of consideration at present the first term, which is supposed 
to represent the diverging waves, and we extend the calculations for E 
(transverse) from the rest of the formula up to C = 3*3 taken at intervals 
of 0*1 for C; we obtain thus the lowest curve given in fig. 11. With the 
help of Table VII, we calculate values of E for depths of about 5, 10, 12, 15, 
and 20 feet, taking in the formula (28) A equal to 

82*5 {1 - 0*14 cos (10*2/c 2 )} 

so that the results apply to the same model at different depths. An example 
of the calculations for one case may be sufficient; Table VIII shows the 
intermediate steps for h = 12*3 ft., L = 20*5. 

Table VIII. 


c 2 . 

c. 

-/3/c. 

E/A. 

e -2-53/c 2 . 

0-68 

0-825 

3-73 

0 *024 

0 -024 

1-12 

1 -06 

2-26 

0*106 

0-106 

1-69 

1-3 

1-5 

0-224 

0-223 

3-25 

1-8 

0-75 

0*508 

0*472 

5-14 

2-27 

0 *374 

0-385 

0-687 

6*8 

2*61 

0 

1 *5 

1 


The results for the five values of h are given in Table IX, and from these 
the curves in fig. 11 have been drawn. 

The general character of the effect of finite depth is clear on inspection of 
the set of curves in fig. 11. If it is required to go to high values of the 
speed-length ratio in a given tank, the ratio of the depth of water to the 
length of the model must be adjusted so that there is no appreciable effect in 
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Table IX. 


h/L. 

x 

c 

0-7 

1-0 

1*5 

1*7 

2*33 

3 


E 

0-7 

7*5 

27 *5 

39 *2 

59 *2 

80 *5 

0-75 

c 

0-7 

1 -o 

1 *2 

1*4 

2 

2*54 

E 

0*7 

7*5 

13 

21 *4 

54 *5 

79 *2 

0-6 

c 

0*7 

1 *0 

1 *3 

1 *8 

2*3 

2*6 

E 

0*7 

7*5 

17 *9 

47 *7 

78 

122 

0*5 

c 

0-7 

1*0 

1 *2 

1 *65 

2*1 

2*38 

E 

0-7 

7*5 

• 13*1 

40 

74 

127 

0-25 . 

c 

0*7 

0*8 

0*84 

1*16 

» 1 *46 

1 *68 

E 

0-7 

1*7 

3*3 

14 *1 

43 *6 

142 *5 



the range of the experiments. Since the curves given here are theoretical 
curves for transverse waves only, each of them ends abruptly at the critical 
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1909.] The Wave-making Resistance of Ships. 299 

Telocity—the resistance being zero after that point. In practice, we know 
that there are no such discontinuities in the resistance curves, and there are 
certain considerations which go to account for this difference. First, as 
regards the transverse waves alone, the preceding formulae show that the 
amplitude tends to become infinite at the critical velocity, although the 
corresponding resistance at uniform velocity remains finite; but, even apart 
from the effects of viscosity, there is a highest possible wave with a velocity 
depending partly upon the amplitude. Secondly, we have left out of 
consideration the diverging waves; but these must become more important 
in the neighbourhood of the critical velocity, for we may regard the two 
systems as coalescing into one solitary wave in the limit as the critical 
velocity is reached. After this point the diverging waves persist, so that 
the effect of these would be of the order of halving the drop in the resistance 
as the critical velocity is passed. 

Finally, we must consider the frictional resistance, which increases steadily 
with the velocity; so that the fall is finally a smaller percentage of the total 
resistance than might appear at first. The curves given in fig. 11 give 
an estimate of a maximum effect of this kind, considering only the transverse 
wave system. 

§7. Further Types of Pressure Distribution . 

The preceding formulae have been built up on the effect of a travelling 
pressure disturbance of simple type; we consider now another type which 
we may use as an illustration. 

Let the pressure system be given by 

p = f{x) = A (h 2 -x 2 )/(x 2 +h 2 ) 2 . 

The type of distribution is graphed in fig. 12. 



Proceeding as in §2, we have 

poo -^2 2 

<f> (k) = 2A I cos kg) dm = ttAkc - ^. 

Jo(® +^) 


(30) 
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300 The Wave-making Resistance of Ships. 

Hence the amplitude of the regular wave-trains formed on deep water in 
the rear of this disturbance is proportional to /c 2 Ae~ Kh , and the effective 
wave-making resistance is proportional to /e 4 A 2 e~ 2Kh . We make the same 
assumption as before, viz., A proportional to v 2 , and write k — g/v 2 \ then the 
resistance is given by 

E = Cv~±e- 2 M V \ (31) 

We use this expression to show how E varies with the constant h of the 
pressure system. Let v = 10 ft./sec., and let E = 1 for h = 0 ; then we find 
the following relative values: 


h. 

R. 

0 

1-0 

1 

0*52 

5 

0*04 

10 

0 *0016 


E decreases very rapidly as h is increased. We have chosen this example 
for the following reason. Consider the motion of a thin infinite cylinder in 
an infinite perfect fluid; if we consider a plane parallel to the direction of 
motion and to the cylinder and at a distance h from it, we find that the 
distribution of excess or defect of pressure due to the motion is of the above 
type. Now, this is not the same as a cylinder moving in deep water at 
a depth h below the free surface, hut it is suggested that as a first approxi¬ 
mation the wave-forming effect is that of an equivalent diffused pressure 
system. The illustration shows how rapidly the wave-making resistance 
diminishes with the amount of diffusion, that is, with the depth h ; this, of 
course, agrees with the experiments on the resistance to motion of submerged 
bodies, and, in fact, with the resistance of submarine vessels. 

In the preceding work no attempt has. been made to connect theoretically 
the constants in the pressure formula with those of the model; since the 
theory rests chiefly on the consideration of transverse waves only, this would 
presumably bring into question the length of entrance, run, and so forth. The 
consideration of any “ transverse ” constants, such as the beam, would need 
a fuller treatment of a diffused pressure system in two dimensions on the 
surface so as to give a more detailed investigation of both transverse and 
diverging wave systems. 
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